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Abstract. In this paper we generalize Jacobi's derivative for- 
mula, considered as an identity for theta functions with character- 
istics and their derivatives, to higher genus/dimension. By apply- 
ing the methods developed in our previous paper [GSMQg] , several 
generalizations to Siegel modular forms are obtained. These gen- 
eralizations are identities satisfied by theta functions with charac- 
teristics and their derivatives at zero. Equating all the coefficients 
of the Fourier expansion of these relations to zero yields non-trivial 
combinatorial identities. 



1. Definitions and notations 

We denote by Tig the Siegel upper half-space — the space of sym- 
metric complex g X g matrices with positive definite imaginary part. 
For e, 5 G (Z/2Z)^, thought of as vectors of zeros and ones, r G Tig 
and 2 G C^, the theta function with characteristic [e, 6] is 



e 



exp m 



S S S S 

(m + |)r(m + |) + 2 *(m + + -) 



A characteristic [e, 6] is called even or odd depending on whether it is 
even or odd as a function of z, which corresponds to the scalar product 
e - 6 G Z/2Z being zero or one, respectively. The number of even (resp. 
odd) theta characteristics is 2^-^29 + 1) (resp. 2^-^29 - 1)). For 
e G (Z/2Z)^ the second order theta function with characteristic e is 



Q[e]{T,z) :=9 



(2r, 2z). 



A theta constant is the evaluation at 2; = of a theta function. 
We drop the argument z = in the notations for theta constants. 
Obviously all odd theta constants vanish identically, and thus there 
are 2^~^(2^-|- 1) non-trivial theta constants with characteristics, and 2^ 
theta constants of the second order. 



First author partially supported by NSF Mathematical Sciences Postdoctoral 
Research Fellowship. 



2 



SAMUEL GRUSHEVSKY, RICCARDO SALVATI MANNI 



A triplet of characteristics [£^2,(52], [£^3,(^3] is called azygetic if 

^_]^^£i-<5i+£2-<52+£3-<53+{£i+e2+e3)-('5i+<52+'53) _ 

A sequence of 2(7 + 2 characteristics [ei,6i], . . . , [£^23+2, ^23+2] is called 
a special fundamental system if the first g characteristics are odd, the 
remaining are even and any triple of characteristics in it is azygetic. 

In the genus 1 case one of the main identities for theta functions is 
Jacobi's derivative formula: 



'11 



ir,z)\ 



2=0 



-710 



{r)e 



There is a long history of possible generalizations of this formula to 
higher genus. We consider g odd characteristics [ei, 61], . 
define their jacobian determinant to be 



[eg,6g], and 



(2) 



D{[ei,6i],...[eg,6, 



71 ^grad 9 



)ir) := 
A grad 9 



A ■ ■ ■ A grad 9 



(r,0). 



Essentially the problem of generalizing Jacobi's derivative formula con- 
sists in expressing some linear combinations of jacobian determinants 
of g distinct odd theta functions as polynomials or rational functions 
in theta constants. When g = 2, such formulas are due to Rosenhain, 
Frobenius, Thomae, Fay, Igusa: we refer to jTh870j . |Fr885j . |Fa79j . 
Ig80| and |SM83j for exact statements and a precise description of the 



situation. We recall from these works that there is a precise conjectural 
formula, which has been proven for g < 5. Moreover, for g < 3 the 
equality 



(3) D{[e^,6^],...,[eg,6g]){T) = ±9 



9+1 



9 



holds if and only if the 2g + 2 characteristics appearing in it form a 
special fundamental system. 

Differential equations for genus 2 theta constants have also been 
studied by Ohyama pi96] and Zudilin [ZuOO]; Grant [H^HH] obt ams a 
nice relation involving only one partial derivative. 

Generalizations of Jacobi's derivative formula in another direction, to 
higher level theta constants in one variable, are derived and discussed 
FKOlj — generalizing these to the higher genus would also be very 



m 



interesting. 

A different generalization of Jacobi's derivative formula involves higher 
order derivatives of theta functions. For example it makes sense in 
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genus 1 to ask for the expression of 
/ 

det 



(4) 



e 



d2 a 



T 



e 



(r,z)|.=o i-^e 
( 



\ 



Am det 



0" 




e 


'l" 


(r) 








'o' 


(r) 


dr 


"l" 









(^,2;) 1^=0 
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as a polynomial in tlieta constants and first-order derivatives with re- 
spect to z of theta functions, evaluated at 2; = (the two determinants 
are equal by the heat equation). We know in fact that 

2 



(5) -e[i](27 



dr 



(e[0](r))/e[l](r)) 



I 

4:71 



dz 



-9 



{r.z)\ 



2=0 



or, equivalent ly, 

e[0](r) 



(6) det d 



47r Idz 



{r,z)l 



To prove this, one can invoke a modular argument, saying that both 
sides are modular of the same weight, thus proportional, and the con- 
stant can be easily computed. Of course such a proof is not very 
revealing, and thus obtaining another proof would be desirable. As 
explained to us by H. Farkas, this identity can also be deduced from 
theorem 5.3 in chapter 2 of [FKOlj by applying Jacobi's triple product 
identity, changing to the argument r/2 and then verifying the resulting 
identity combinatorially for each coefficient of the Fourier series. 

In this paper we shall generalize both the identities (0) and © to 
higher genus. It would be interesting to understand the combinatorial 
meaning of these generalizations similarly to the one- variable identities 
above or obtain alternative combinatorial proofs, but these questions 
lie beyond the scope of the current work. 

The main tool will be a consequence of Riemann's addition theorem 
relating the first ^-derivatives of odd theta functions with characteris- 
tics to the second 2;-derivatives of second order theta functions. This 
has also been the main tool in our paper |GSM03j . where we showed 
that generically a principally polarized abelian variety is uniquely de- 
termined by the gradients of odd theta functions at 2; = 0. 

Remark 1. We note that the classical generalization of Jacobi's de- 
rivative formula can be given an interpretation in terms of theta series 
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with harmonic polynomial coefficients. In fact monomials of degree 
(7 + 2 in the theta constants are theta series relative to the quadratic 
form 4- lg_|_2 and harmonic polynomial "1" . The jacobian determinants, 
on the other hand, are theta series relative to the quadratic form 4 • 1„ 



and harmonic polynomial "det", Ig83 



This is the simplest pair of theta series with harmonic polynomial 
coefficients. Our generalizations can also be interpreted in this way. 
For example in genus one, while, as we wrote, in the ffist two cases the 
harmonic polynomials are 1 and x, in our case it is a polynomial in two 
variables: — y"^- 

As a further consequence of our formulas, we shall give a characteri- 
zation of the locus of reducible principally polarized abelian variety in 
terms of vanishing of certain derivatives of odd theta functions. 



2. The symplectic group action 
Let Tg := Sp{2g, Z) be the integral symplectic group; it acts on TCg 



by 



where M 



M ■t:={At + B){Ct + D)-^, 



A B 
C D 



G r„. A period matrix r is called reducible if 



there exists M G such that 



M-r 



Tl 
T2 



, Ti G Hg^, 91+92 = 9] 



otherwise we say that r is irreducible. 

We define the level subgroups of the symplectic group to be 



Tg{n) :-- 



M 



Vg{n, 2n) := {M G Vg{n) \ diag(^*5) = diag(C*D) = mod 2n] . 

A function F : Tig ^ C is called a modular form of weight k with 
respect to T G Tg if 

F{M ■ t) = det(Cr + DfF^r), 'iM G F, Vr G Ug 
The theta functions transform under the action of F^ as follows: 



M 



{M-T, \Ct + D)-'z) 



{e, S, M, r, z)det{CT + D)^9 
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where 



M 



D -C 
-B A 



£\ ldmg{C'D) 
+ ' diag(A*fi) 



taken modulo 2, and 5, M, r, z) is some complicated explicit func- 
tion. For more details, we refer to |Ig72| and |R,F74j . 

Theta constants with characteristics are modular forms of weight 1/2 
with respect to rg(4, 8). In this case (5, M) := 6, M, r, 0) is an 
eighth root of unity that does not depend on r. 

Differentiating the theta transformation law above with respect to 
some Zi and then evaluating at z = 0, we see that 



dz, 



■e 



M 



d 



(e, 5, M) det(Cr + Df''' ^^(Cr + D),j—e 



{r,z)\ 



z=0- 



Denoting by grad 9 



spect to ^1,^2, 



^ Zg ^ J Z 



(r) the gradient of the theta function with re- 



0, the above formula becomes 



grad 6 



M 



(M-r) =0(e, 5, M) dei{CT+D)^{CT+D)grad 9 



ir). 



As a consequence, the jacobian determinant D{[ei, 6i], . . . [sg, Sg]){T) 
is a modular form of weight ^g + I with respect to Tg{4:, 8) (see |Ig80| 
and jSMEHI). 



3. Some multilinear algebra 

For our purposes we need some results from linear algebra, which 
we recall and prove for the sake of completeness. We are grateful to 
C. De Concini, A. Maffei, D. Zagier and one of the referees for useful 
suggestions about these topics. 

To any A G MatgxgiC) we associate the {g — 1) x [g — 1) matrix 
A whose entries are the determinants of 2 x 2 minors of A obtained 
taking the first line and the first column and letting the other row and 
column vary, i.e. 

Ai, ■= det Ali , 2<i,j <g. 

We observe that decomposing the matrix A in blocks 

/an 'z\ 
^-[w B)^ 
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A 



with 5 a ((7 — 1) X ((7 — 1) matrix, and z,w E'C^ ^, we have 

Aij := det(aiii? — w^z) 
With these notations we have 
Lemma 2. 

(7) a{\^ det A = detA 

Proof. It is trivial when Oi 1 = and it is an immediate consequence of 

1 \ fail 

when ai 1 7^ □ 

We denote := g{g + l)/2 and for any G O we denote by its 
symmetric tensor square. Then the following is true 

Lemma 3. Let vi, . . . ,vn G C^- Then 

{N\)vl AvlA---Avlf= sign(s)(f,(i) A t;,(2) A ■ ■ ■ A v,(g))- 

seSM 

{Vs{l)AVs{g+l)AVs{g+2)^- ■ ■^Vs{2g-l))-{Vs{2)^Vs{g+l)AVs{2g)^- " "Af 5(39-8) ) 
{Vs(3)^Vs{g+2)AVs(2g)^- " ■Afs(4g-6)) • • • {Vs{g)^Vs(2g-l)^Vs{3g-3)^- " -^VsiN)) 

Proof. Since the LHS is SL{g,C) invariant, it can be expressed as a 
polynomial in determinants of g x g minors of the g x N matrix with 
columns being f 's. Moreover, this polynomial must be homogeneous 
of degree g + 1 in these determinants, each Vi has to appear in it ex- 
actly twice, and it has to be skew-symmetric in v's. For this reason 
first we sum over all possible permutations with the signs. We further 
observe that if the same two vectors appear in two different determi- 
nants, then the expression vanishes. Thus the expression has to be 
a sum of monomials each of degree g + 1 in the determinants, such 
that each v appears exactly twice, and no pair of f 's appears twice in 
two different determinants. Thus the expression is forced to be exactly 
that of the statement, up to a multiplicative constant, which is easily 
computed. □ 

4. 9's AND 9'S 

A special case of Riemann's bilinear addition theorem for theta func- 
tions (see |Ig72| , |RF74j . |Mu84j ) is 



(8) e[a](r,z)e[a + e](r,0) = l Yl ("l)""^^ 



29 



{r,z)9 
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which is vahd for all r and z. Taking a sum of these with different 

signs, we get, for any 5 G {7LlT£)^ 

(9) 

E (-ir^eH(r,^)0[a + e](r,O) = 

«e(z/2Z)9 



- r 



We assume that the characteristic [e, 5] is odd, differentiate this re- 
lation with respect to Zi and Zj^ and then evaluate at z = 0. Denoting 
by Qiehij) the g ^ g symmetric matrix with entries 



(r,0)ai^ 



(r,0). 



and by K^^ir) — the g ^ g symmetric matrix with entries 

A.,, .,(r) := (9,,9,^e[5](r,0))e[£](r,0) - (9.,9,^,e[£](r, 0))e[5](r, 0), 
we thus have (see |GSM03j ) — notice that C^^ = if [e, 5\ is even 
Lemma 4. 



(10) 



5^ 



and the "inverse" 
Lemma 5. 
(11) A 



— y 

29-1 



-l)-^a,(r) 



/3g(Z/2Z)9 



We remark also that 
(12) Q5(r) = 2^rarf ^ 



(r) ^gradO 



[T . 



5. Generalized Jacobi's derivative formulas 

To generalize the first result of the introduction, we introduce the 
matrix-valued differential operator 



I _d_ 1 d 

1 d a 

2 dT21 dT22 



\ 1^ 



1^ \ 

2 drig » 

1 9 

2 dT23 

drgg ) 



Then we have 
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Theorem 6 (First generalization). For any e ^ 5 the following holds: 

ce[5f3 det(p(e[e]/e[5]) 

(13) = E (-l)^-(-n+-+->.)z}([e + 5, a,J, ...[e + 5, a,^]f 

{ttij |[e+(5,Oi^] odd} 

for some computable constant c. 

Proof. For any characteristics e, 5 we have by definition 
Q[5fV{Q[e]/Q[5]){r) = AmAe5{T). 

Thus 

e[5]'3det {V{Q[e]/Q[5])){T) = {Aniy det{Ae5)ir). 
Now, using the resuh of lemma 2, we get 



det Assir 



det 



1 



29- 



-1 



E 



det 



29 



-1 



E 



{a|[e+(5, a] odd} 



e + 5 
a 



(r)* gradO 




{Q:|[e+<5, a] odd} 

When we expand this determinant, each summand will be of the type 



sign(/i) 



sign(cr) 9,^(1)6' 



£ + 5 

ai 

'e + 5 
ai 



6 + 6 



an 



(r) 



e + 6 
an 



(r) 



for some permutations a and fi. Taking the sum of these for all pos- 
sible permutations a and /i gives exactly the square of the jacobian 
determinant, so that we end up with 

proving the theorem. □ 

At this point, we observe that our relations are not trivial. In fact 
each term appearing in the RHS is not identically zero, jSM83j . We 
remark that the set of characteristics appearing in the jacobian de- 
terminant above is syzygetic, while in all the other generalizations of 
Jacobi's derivative formula only azygetic sets appear, cf. |Ig80| . 

If we would like to have relations involving the derivatives of the 
second order theta constants with respect to Tij, then since the matrix 
C^si'^) has rank 1, we have the following 
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Proposition 7. If g > 2, then 



(14) 




\a-S 



A. 



0. 



For the LHS of (fT!^ . we have the following non- vanishing: 
Theorem 8. For all possible pairs e ^ 6 & (Z/2Z)^ the expression 

e[d]^'det {v{e[e]/e[6])) 

is not identically zero in t. 

Proof. We shall prove a slightly more general result: that for any pair 
of distinct even characteristics [e, a] and [6, /?] the expression 

29 



e 



det [Vie 



is not identically zero. 

Indeed, we know that the symplectic group Tg acts doubly tran- 
sitively on the set of even characteristics, and we have the following 
transformation formula 



V e 



M 



a 



M 



{e,a,S,p,M){CT + Dyv 9 



{M-t) 

t){Ct + D) 



with (f){e, a, 5, /5, M) an eighth root of unity. 

We learnt from JBZOSJ that for some specific [e, 0] and [5, 0] we have 



not identically zero; since all such expressions are 



permuted by the symplectic group action, they are all not identically 
zero. □ 



For the second generalization of Jacobi's derivative formula, for any 
set of + 1 := \g{g + 1) + 1 characteristics Eq, Si, 
the matrix 



M^Eq.Ei, . . . ,£:Ar)(r) 



, we introduce 
Q[eN]\ 



d 
drac 



'0[eN]J 



10 



SAMUEL GRUSHEVSKY, RICCARDO SALVATI MANNI 



It is a well-known fact (see |Sa83j ) that det M{eo,ei, . . . ,eN){T) is a 
modular form of weight ^{g + 2){g + 3) relatively to Tg{2, 4). Moreover, 
if it is not identically zero, then the theta constants 

e[£o],e[£i],...,e[£jv] 

are algebraically independent. We set 

'■= £0 + £k 

and let Cepir) be the vector in with entries Ce^^ij{r) (before we 
thought of C as a matrix, but now we write down all the matrix ele- 
ments in a single vector). Using the results of Lemmata 4, 2 and 1, we 
get the following 

Proposition 9 (Second generalization). For some computable constant 
c 

ce[eof-Uet M{eo,ei, . . . ,eN){r) 

{-iy''^'^det{Cs,/3,ACs,p,A---ACs^p^) 

/3i,...,/3jvg{Z/2Z)9 

By lemma 2 the RHS can he expressed as a homogeneous polynomial of 
degree g + 1 in jacobian determinants. 

Remark 10. Recalling the definition of theta functions, we can ex- 
pand all of the above identities in Fourier series in r or equivalently 
in power series in qij := exprjj. The coefficients of these expansions 
will then be some rather complicated but quite natural combinatorial 
quantities in several variables, and the equality of the RHS and LHS of 
any of the above would then yield a non-trivial multidimensional com- 
binatorial identity, which it would be interesting to understand and 
prove combinatorially. 

6. An application in genus 2 

We will now work out in detail the situation in the case of genus 2. 
Indeed let us write down (fT^ for e = [00] and 6 = [10]: 

(15) e[io]^ det(D(e[oo]/e[io]) = cD([io, lo], [lo, ii])^, 

with c a known constant. 

Using lemma 0] to express the RHS in terms of theta constants of 
the second order and their derivatives, we get (we denote dij := d^^ to 
simplify notations) 

(6 [00] an e [10] - 0[io]aiie[oo]) (e[oo]a22e[io] - e[io]a22e[oo]) 

- (e[oo]ai2e[io] - e[io]ai2e[oo])' 
+ (e[oi]ane[ii] - e[ii]ane[oi]) (e[oi]a22e[ii] - e[ii]a22e[oi]) 
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Clearly we get the same equation if we chose e = [01] and 6 = [11]. 
Thus in all we get three different equations. 

In |BZn3j it is shown that in genus 2 there are 2 ■ 2^ + 2 = 10 
algebraically independent quantities among 4 theta constants of the 
second order and their 4 ■ 3 = 12 derivatives. Thus there are 6 non- 
trivial algebraic relations among theta constants and their first-order 
derivatives. So far we have obtained three such equations, and three 
more can be obtained by writing down formula (fT^ . 

In |BZ08j some other 6 independent relations are given. We shall 
prove that they are all consequences of (fT^ and Q. 

Indeed, for genus two formula Q reads 



Applying Q to rewrite the RHS in terms of theta constants of the 
second order, we finally see that 



is a polynomial in the theta constants of the second order. This equa- 
tion is up to a rational function equal to one of the equations in |BZ03j . 
From the other choices of characteristics e and 6 we get the other 5 
equations. We observe that these 5 equations can also be obtained 
from the first one by the action of In this spirit we remark that in 
genus 1, we have four variables and the relation is an immediate conse- 
quence of (IH) and (0). These relations can also be obtained considering 
the determinant of (fT^ . 

This method allows us to give a conjectural description of the sit- 
uation in the genus 3 case. By the results of [BZ03| we know that 
among 56 variables (all Q[e] and their derivatives 5jj0[£]) there are 21 
algebraically independent ones. Thus there are 35 algebraic relations. 
We know that there is a unique polynomial relation among the Q[e], 



of degree 16, cf. |vGvG86] : let us denote it by R{t). Thus we have 

R{t) = 0, dijR{r) = 



and other 28 relations obtained by applying (jl4|) . since in genus 3 there 
are exactly 28 odd characteristics. So in total we have 35 relations that 
we conjecture to be algebraically independent. 

7. Characterization of the reducible locus 

We finish the paper by giving a characterization of the locus of re- 
ducible abelian varieties. Different characterizations of the reducible 



(e[oi]9i2e[ii] - e[ii]9i2e[oi])' = o 




e[io]^det(p(e[oo]/e[io]) 
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locus are known: in |EL97j it is characterized in terms of the dimen- 
sion of the singular locus of the theta divisor, in |Sa83j — in terms of the 
non-maximality of the rank of matrix P(t) with rows -^Q[€]), 



and columns corresponding to all e G (Z/2Z)^, in IMII — in terms 
of the vanishing of certain theta constants. Here we use the vanishing 
of certain first derivatives of theta functions evaluated at zero. 

Proposition 11. A ppav with a period matrix t is reducible if and only 
if there exist some M e and some k < g such that if we write any 
odd characteristic [e, 6] as [ei 62, Si 62], where [ei, 61] is a k-dimensional 
characteristic, and [£2,^2] "is {g — k) -dimensional, then 







for all i < k for [ei, 61] even, and for all i > k for [ei, 61] odd. 



Proof. Suppose the period matrix M-t splits as M-r 



with 



n 

T2^ 

Ti G Tik and T2 G 'Hg-ki so that the theta functions with characteristics 
factor as follows: 



e 



{M-T,z) = e 



^1,^1, 



■6 



The vanishing of the derivatives in question is immediate for M ■ t 
by differentiating and evaluating at 2; = 0; thus the "only if" part is 
proven. For the "if" part, assume the vanishing of derivatives as stated. 
Then according to (fTTj) we have Ces,ij = and consequently A^^^j^ = 
for l<i<k<j<g and all odd [e, S]. Thus the matrix P{r) does not 
have maximal rank, and thus, by the results of |Sa83j . it corresponds 
to a reducible abelian variety. □ 
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